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PREFACE 


This  project  was  conducted  in  support  of  the  Armstrong 
Laboratory  Human  Resources  Directorate,  Aircrew  Training  Research 
Division,  research  concerning  simulation  of  infrared  imaging.  This 
effort  was  managed  by  the  University  of  Dayton  Research  Institute 
(UDRI) ,  under  Contract  F33615-90-C-0005 ,  Work  Unit  2743-25-17 
Flying  Training  Research  Support.  The  laboratory  contract  monitor 
was  Patricia  A.  Spears. 

This  effort  supports  the  training  technology  objectives 
described  in  the  present  Armstrong  Laboratory  Research  and 
Technology  Plan.  The  goal  of  the  present  research  was  to  provide 
an  understanding  of  the  computational  aspects  of  the  infrared 
imaging  simulation  and  clarify  the  shortcomings  in  the 
computational  model  presently  being  used. 

I  wish  to  express  my  thanks  to  Ms.  Marge  Keslin  for  final  edit 
of  the  manuscript. 


SIMULATION  IN  INFRARED  IMAGING:  USING  ELECTRICAL  CIRCUIT 
PRINCIPLES  TO  MODEL  HEAT  TRANSFER 

SUMMARY 

One  approach  to  computer  simulation  of  infrared  (IR)  imaging 
is  to  model  heat  transfer  between  media  using  electrical  circuit 
theory  principles.  A  previous  implementation  of  such  an  approach 
produced  significant  errors  in  surface  temperature  predictions  for 
database  features.  The  purpose  of  this  effort  was  to  reveal  the 
source  of  these  errors  and  explore  alternatives. 

Heat  transfer  is  continuous  across  both  space  and  time.  For 
individual  objects  or  features,  heat  transfer  and  storage  are 
highly  dependent  upon  individual  material  properties,  exposure  of 
features  to  the  sun's  radiation,  and  thermal  insulation  from 
objects  beneath  or  adjacent  to  the  feature  of  interest.  Such 
characteristics  make  the  building  of  general  models  quite 
susceptible  to  significant  errors. 

Analyses  from  this  report  revealed  that:  (a)  some  temperature 
prediction  errors  could  be  attributed  to  the  discrete 
representation  of  database  features  and  failure  of  the  model  to 
characterize  thermal  insulation  for  individual  database  features, 
and  (b)  inclusion  of  capacitance  or  heat  storage  capabilities  in 
the  model  was  problematic  because  material  coefficients  specific  to 
individual  features  led  to  inconsistencies  in  the  prediction  of 
daily  heating  and  cooling  cycles  for  database  features. 

No  simple  solution  to  these  problems  exists.  The  desire  to 
simulate  the  infrared  signature  of  a  large  number  of  database 
features  under  a  wide  range  of  conditions  necessarily  requires  a 
model  of  great  complexity.  The  results  of  the  present  analyses 
provide  a  focus  for  future  improvements  in  the  model. 

INTRODUCTION 

Forward-looking  infrared  (FLIR)  systems  used  on  military 
aircraft  are  an  integral  part  of  the  weapon  systems  package. 
Infrared  imagery  from  the  FLIR  sensor  is  typically  displayed  on 
cathode-ray  tubes  (CRTs) ,  appearing  somewhat  similar  to  black  and 
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white  television.  To  simulate  the  display  of  the  FLIR  image  on  the 
CRT,  the  scene  irradiance  from  the  database  image  impinging  on  the 
FLIR  sensor  is  estimated,  and  computational  procedures  are  employed 
to  translate  estimated  irradiance  into  computer  digital-to-analog 
code  (DAC)  values  for  display  on  a  CRT. 

Previous  analysis  of  the  IR  imaging  simulation  software  (Evans 
&  Crane,  1990)  in  use  at  the  Armstrong  Laboratory,  Human  Resources 
Directorate,  Aircrew  Training  Research  Division  (AL/HRA) ,  at 
Williams  AFB,  AZ,  documented  inadequacies  in  the  software 
predictions.  Specific  examples  of  shortcomings  were  the 
variability  in  diurnal  feature  surface  temperatures  and  the 
inability  of  database  objects  to  store  or  dissipate  heat  (i.e., 
thermal  capacitance  of  objects)  over  time. 

Figure  1  is  a  sample  graph  taken  from  Evans  and  Crane  (1990) . 
This  graph  denotes  surface  temperatures  for  5  sample  database 
features  as  a  function  of  time  of  day.  Interestingly,  this  graph 
is  the  wide  fluctuation  in  the  surface  temperature  of  soil  as  the 
day  progresses.  For  example,  soil  is  much  cooler  than  the  other 
features  during  the  early  morning  hours  but  heats  up  much  more 
relative  to  the  other  features  during  the  day.  The  FLIR 
impression,  then,  of  the  soil  against  a  metal-based  object  during 
midday  would  be  that  the  soil  is  hotter  than  the  metal  object. 
Given  a  scenario  where  pilots  are  using  FLIR  to  locate  or  detect  a 
metal  structure  during  the  day,  their  objective  would  be  to  find  a 
hot  area  (i.e.,  the  metal)  in  the  midst  of  a  cooler  background 
(i.e.,  the  soil).  Therefore,  such  erroneous  predictions  by  the 
software  are  of  considerable  concern. 


REPRESENTING  HEAT  TRANSFER  THROUGH  ELECTRICAL  CIRCUITS 

The  distribution  of  thermal  energy  in  a  system  may  be 
represented  through  analogy  to  an  electrical  circuit.  Sources 
(e.g.  ,  2the  sun  or  sky)  of  radiant  energy  or  flux  measured  in 
watts/m  may  be  represented  as  current  within  ^n  electrical 
circuit.  Conductance  (or  alternatively  conductance  1  =  resistance) 
of  heat  across  media,  measured  in  watts/m /C  (or  °C/watt/m2  for 
resistance) ,  may  be  represented  as  conductance  (resistance)  in  an 
electrical  circuit,  and  temperature  differences  in  a  thermal  system 
may  be  represented  as  voltage  differences  across  elements  of  a 
circuit  (Suits,  1985) . 


Figure  2  denotes  one  circuit  representation  of  the  steady- 
state  temperature  for  a  database  object  in  thermal  equilibrium. 
Input  flux  or  current  sources  consist  of  the  sun  (PS1^)  and  the  sky 
(Ptky)  .  Radiation  emanating  from  the  database  object  of  interest, 
^object *  is  estimated  through  the  Stefan-Boltzmann  equation: 


^object 


(1) 


2 


1  7  13  1* 

4  10  16  22 
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Figure  1.  Surface  Temperature  as  a  Function  of  Time  of  Day 
(Arizona  in  August) . 


Steady-State  Feature  Temperature  Circuit. 
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where  e  is  average  material  emissivity,  a  =  5.6686  X  10  W/ (m  K  ), 
and  T  is  temperature  in  °Kelvin.  This  4th  order  equation  of 
temperature  is  approximated  in  a  linear  fashion  for  circuit 
representation  such  that: 


P  =  eaT4  =  e(316+  4.62T)  (2) 

where  T  in  the  linear  approximation  is  now  in  degrees  Celsius.  The 
circuit  in  Figure  2  shows  2  branches  representing  the  2  flux 
components,  316e  and  4.62eT  (or  4.62cV  where  voltage  is  substituted 
for  temperature) ,  for  the  emission  from  the  database  object  itself. 
Other  branches  of  the  circuit  in  Figure  2  represent  independent 
temperature  sources  (rain,  air,  support  base) .  Each  source  is 
modified  by  a  conductance  (resistance)  across  the  other  branches 
(i.e.,  hk,  hc,  and  Rc;  the  conductance  across  the  feature  of 
interest,  from  the  wind,  and  from  rainfall) .  Finally,  in  the 
steady-state  or  time  invariant  approach,  the  temperature  (voltage) 
across  the  database  object  of  interest  is  equivalent  to  the 
resistance  of  the  database  object  multiplied  by  the  flux  through 
the  database  object  and  added  to  the  temperature  (voltage)  of  the 
base  support  (T^J  . 

Fourier's  Law  describes  the  rate  of  heat  flux,  P,  through  an 
object  by  the  formula: 

P  =  (k/d)  (T1  -  T2)  in  watts/M2  (3) 

where  k  is  the  thermal  conductivity  of  the  material  in  watts/ 

(M  °K)  ,  d  is  the  one-dimensional  depth  or  thickness  of  the 
material,  and  T,  and  T2  denote  the  temperature  in  °K  on  either  side 
of  the  material.  P  is  the  flux  which  is  analogous  to  current  in 
our  electrical  circuit,  the  temperature  difference  (T,  -  T2)  in 
Fourier's  Law  is  analogous  to  the  voltage  difference  across  a 
medium  in  an  electrical  circuit,  and  (k/d)  represents  conductance 
in  either  situation.  Note  that  Figure  2  contains  the  resistance 
d/(2k)  and  not  d/k  as  shown  in  Fourier's  Law.  This  change 
represents  the  amount  of  resistance  across  one-half  the  depth  of 
the  material  and  enables  us  to  remain  consistent  with  further 
developments  of  the  model  in  this  report  when  T^^  will  be  replaced 
with  a  capacitor.  The  conductance  parameters  provide  the  rate  of 
exchange  per  unit  of  temperature  difference  across  objects 
contained  in  our  closed  system  (e.g.,  database  object,  air,  base 
support,  rain) . 

The  circuit  in  Figure  2  is  a  thermal  model  of  the  world  from 
the  perspective  of  tha  database  feature  of  interest.  Actual 
temperature  gradients  are  continuous  and  all  materials  have  some 
ability  to  capture  a  id  distribute  thermal  energy.  Figure  2,  then, 
is  necessarily  a  simplification  of  the  world  from  the  database 
object's  point  of  view.  Along  with  this  simplification,  certain 
quirks  of  the  representation  can  play  an  important  role  in 
determining  surface  temperature  predictions.  For  example,  when 
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there  is  no  temperature  (voltage)  or  radiant  flux  (current)  input 
from  any  of  the  sources,  the  temperature  or  voltage  difference 
across  the  branches  is  zero.  Zero  voltage  denotes  a  surface 
temperature  of  0  °C.  Iherefore,  in  the  current  implementation  of 
the  model,  surface  temperatures  migrate  toward  freezing  when  the 
current  and  voltage  sources  contained  in  the  model  are  set  to  zero. 
Feature  surface  temperature  migration  toward  freezing  is  an 
important  artifact  of  the  model  and  shall  be  discussed  in  more 
detail  later  in  the  report. 

The  solution  for  feature  surface  temperatures  in  the  original 
software  is  accomplished  by  using  nodal  and  mesh  analysis  from 
circuit  theory  to  solve  for  Tss  or  voltage  in  Figure  2.  The 
solution  of  Tss  in  Figure  2  is: 

Psun  +  Psky  -  3166  +  h£Tajr  +  h.T^  +  RcTrajn 

Tss  - - Celsius  (4) 

4 . 6  e  +  hc  +  hk  +  Rc 

hc,  hk,  and  Rc  represent  conductances  for  wind,  the  database  object 
of  interest,  and  rain,  respectively.  Tajr,  T^^,  and  Train  represent 
temperature  sources  for  air,  a  base  of  support,  and  the  rain,  all 
measured  in  °Celsius.  Further  explanation  of  Equation  (4)  can  be 
found  in  Evans  and  Crane  (1990)  while  the  general  approach  can  be 
found  in  Suits  (1985) .  Equation  (4)  denotes  a  steady-state 
solution  for  feature  surface  temperatures,  and  Figure  1  denotes  a 
plot  of  Equation  (4)  for  the  sample  features.  The  steady-state 
solution  applies  when  database  features  do  not  store  or  build  up 
heat  capacity  for  release  at  a  later  time.  The  calculation  of 
surface  temperature  above  serves  for  both  small  and  massive  objects 
in  the  current  version  of  the  software  used  at  AL/HRA. 

Before  proceeding  with  modifications  of  Figure  2,  it  proves 
useful  to  examine  a  common  scenario  for  the  circuit  in  Figure  2 
which  will  highlight  problems  inherent  in  the  model.  For  a  calm 
day  with  no  wind  or  rain,  hc  and  Rc,  the  air  conductance  and  the 
rain  conductance,  will  be  zero.  Zero  conductance  amounts  to 
infinite  resistance,  allowing  no  current  to  flow  through  the 
respective  branches  in  Figure  2.  The  current  or  flux  generated  by 
the  three  independent  generators,  Pst^  +  Psky  -  316e  is  then  divided 
across  2  branches,  the  first  with  resistance  1/(4.626)  and  the 
second  with  resistance  d/(2k)  and  independent  temperature  (voltage) 
source  Tbase  which  is  generating  current  itself.  The  solution  for 
Tss  is  then: 


Tss 


Psun  +  Psky  "  3166  +  h.T^ 


4.626  +  hi. 


Celsius  (5) 


psun  and  psky  represent  radiant  flux  from  the  sun  and  sky  modified  by 
the  absorption  coefficient  of  the  specific  database  feature. 
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Radiant  flux  from  the  sky  is  estimated  using  the  Idso-Jackson 
Equation  (Suits,  1985) : 

Psky  “  aTair4Cl  "  *261  exp{-7.77  X  10'4  X  (273-Tair) 2}  ]  (6) 

o  ”82 

wjiere  Tajr  is  air  temperature  in  Kelvin  and  a  =  5.6686  X  10’  w/ (m 
K  )  is  the  Stefan-Boltzmann  Constant.  Figure  3  shows  a  plot  of  P8ky 
as  a  function  of  air  temperature  in  °F.  A  typical  value  for  Psky, 
not  attenuated  by  material  absorption,  is  300  watts/ni 
(60  °F)  . 

The  absorption  coefficients  for  each  material  vary  as  a 
function  of  wavelength,  making  it  difficult  to  represent  absorption 
for  each  material  with  a  single  scalar.  However,  for  the  sample 
features,  soil  typically  had  the  lowest  absorption  with  an  average 
of  .4  and  water's  aosorption  coefficient  was  typically  about  .99. 
Figure  4  shows  irradiance  absorbed  from  the  sun  and  sky  under 
specific  conditions  (Arizona  in  August,  no  wind  or  rain)  for  each 
of  the  sample  features.  While  the  contribution  of  the  sky  was 
estimated  using  the  Idso-Jackson  Equation,  the  contribution  from 
the  sun  was  estimated  using  tables  from  Suits  (1985) . 

In  Figure  4,  irradiance  absorbed  by  the  sample  features  varies 
from  approximately  85  to  145  watts/mz  at  night  to  a  maximum  between 
590  and  980  during  midday.  Psun  +  Psky  -  316c  (c  being  between  .9 
and  1  for  all  features)  could  be  as  low  as  -215  for  soil  at 
nightti  ae  or  as  high  as  680  for  water  during  midday.  From  this 
range  of  values  in  Equation  (5)  ,  we  add  the  product  h^^^.  The 
conductance,  hk,  across  an  object  to  a  base  support  with 
temperature,  T^^,  ranges  from  a  minimum  of  2.7  for  soil  to  a 
maximum  of  1052.6  for  the  bridge  crane  in  the  sample  features.  For 
soil,  then,  the  numerator  in  Equation  (5)  is  approximately  -215  + 
2*7Tbase  at  where  T^^  is  given  in  degrees  Celsius.  For  Tbase 

less  than  80  °C  (176  °F)  ,  then,  the  soil  temperature  will  be  less 
than  0  °C  or  freezing  temperature. 

Figure  1  shows  the  result  of  this  effect  where  soil  surface 
temperature  dips  quite  low  at  night.  In  addition,  because  hk  =  2.7 
for  soil,  which  is  quite  low  for  our  sample  features,  increases  in 
Psun  during  the  day  in  Equation  (5)  will  have  a  more  dramatic  effect 
on  the  change  in  soil  surface  temperature  than  the  change  in  the 
surface  temperature  of  other  features.  This  effect  is  also 
emphasized  in  Figure  1. 

Given  these  facts,  we  find  that  Tss  wil^  be  heavily  dependent 
on  hk,  the  conductance  (measured  in  watts  M  1  °C)  across  a  feature 
(hk  --  2k/d  being  halfway  across  a  feature  and  k/d  across  the  entire 
feature)  which  denotes  the  passage  of  flux  to  a  base  of  support  to 
which  the  feature  is  attached.  As  hk  increases,  Tss  approaches  T^,^. 
For  the  5  sample  features  in  Figure  1  (bridge  crane,  aircraft 
parking,  soil,  water,  and  evergreen  tree),  hk  was  1052.6,  7.9,  2.7, 
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Figure  4.  Sun  and  Sky  Irradiance  Contribution  to 
Feature  Surface  Temperature  Model. 


19.1,  and  25  (w  M*1  °C) ,  respectively.  For  the  bridge  crane  where 
hk=1052.6,  TS5  will  always  be  approximately  T^^.  This  version  can 
be  seen  by  inserting  sample  values  into  Equation  (5)  .  In  the 
original  version  of  the  software,  T^^  was  set  to  the  air 
temperature  and  Evans  and  Crane  (1990)  found  that  the  surface 
temperature  of  the  bridge  crane  always  mimicked  the  current  air 
temperature  quite  closely. 

The  example  above  points  out  a  shortcoming  in  the  model  from 
Figure  2.  Metal  in  direct  sunlight  is  able  to  attain  temperatures 
much  greater  than  the  surrounding  air  temperature.  In  the  model, 
however,  the  surface  temperature  of  materials  which  are  highly 
conductive  (e.g.,  metals,  aluminum)  will  approach  T^^  when  the 
circuit  from  Figure  2  or  Equation  (4)  is  employed.  Valid  thermal 
modeling  of  these  materials  requires  that  they  be  allowed  to  store 
energy  in  some  fashion. 

Alternatively,  soil  is  a  very  bad  conductor  with  a  conductance 
of  hk  =  2.7.  Figure  1  shows  large  variability  throughout  the 
diurnal  temperature  cycle  of  soil  as  mentioned  in  the  last 
paragraph.  For  soil  as  represented  in  the  circuit  of  Figure  2, 
then,  with  very  little  conductance  to  a  base  support,  there  is  very 
little  external  input  to  maintain  a  reasonable  temperature  (i.e., 
anything  above  freezing  or  zero  voltage)  when  flux  input  from  the 
sun  is  near  zero.  In  this  specific  instance,  we  realize  that 
ground  voltage  for  soil  (or  thermal  support  from  beneath  the 
surface)  should  probably  be  set  near  60  °F.  A  simple  solution  for 
soil  would  be  to  set  the  ground  voltage  to  +60  (15.6  °C)  which 
would  yield  a  circuit  solution  of: 


TSS  = 


^sun  +  ^sky  ~  3166  +  hkTb<le 

-  +  15.6  Celsius 


4.6c  +  hb 


(7) 


which  simply  adds  15.6  °C  (60  °F)  to  the  steady-state  temperature. 
This  solution  will  not  suffice,  though,  as  a  general  solution 
because  (a)  it  is  not  desirable  to  add  60  °F  to  all  feature  surface 
temperatures,  and  (b)  the  same  diurnal  variation  in  soil  surface 
temperature  would  exist  which,  as  shown  in  Figure  1,  is  quite  high. 


Alternative  solutions  to  this  problem  involve  either  (a)  the 
addition  of  other  temperature  (voltage)  sources  which  conduct  heat 
to  and  from  our  feature  of  interest  or,  (b)  allow  our  feature  of 
interest  to  store  and  release  heat  over  time.  Alternative  (b)  has 
been  studied  more  closely  to  date,  so  we  will  begin  further 
analysis  from  this  point  of  view  and  include  modifications  of 
alternative  (a)  along  the  way. 
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THERMAL  CAPACITANCE  MODELING 


To  allow  database  objects  to  store  and  release  thermal  energy 
over  time,  the  circuit  analogy  must  include  a  capacitor.  A  simple 
circuit  representing  database  objects  with  heat  storage  and 
releasing  capabilities  was  suggested  by  Hirsch  (Personal 
Communication)  similar  to  Suits  (1985).  Figure  5a  is  a  duplicate 
of  Figure  2  with  the  exception  that  the  temperature  source  Tbase  has 
been  changed  into  a  capacitor  which  represents  the  temperature 
storage  capability  of  the  database  object  itself.  Figure  5a  is  a 
circuit  representation  for  massive  database  features.  It  does  not 
include  a  base  of  support. 

Using  the  principles  of  linearity  and  superposition,  the 
solution  for  feature  surface  temperature,  T,  in  Figure  5a  can  be 
broken  into  two  components,  a  steady-state  or  forced  component  and 
a  time-varying  or  natural  component.  By  replacing  the  capacitor  in 
Figure  5a  with  an  open  circuit,  the  steady-state  temperature  may  be 
determined.  Using  circuit  theory  principles,  steady-state  surface 
temperatures  for  massive  bodies  (bodies  which  are  at  least  a  few 
feet  thick)  can  be  expressed  as: 

^sun  +  ^sky  “  3166  +  ^c^air  +  ^c^rain  0 

T  = -  Celsius  (8) 

4.6c  +  hc  +  Rc 

The  term  "steady-state"  used  here  is  somewhat  of  a  misnomer.  The 
flux  input  from  the  sun  and  sky  vary  in  a  sinusoidal  fashion 
throughout  the  day  and,  therefore,  Tss  will  not  remain  constant. 
Tss  is  the  result  of  applying  the  forcing  function  (i.e.,  the  flux 
input)  to  the  circuit  in  Figure  5a  when  no  current  flows  down  the 
path  of  the  capacitor. 

Since  no  flux  flows  down  the  path  of  the  capacitor,  Tss  in 
Equation  (8)  does  not  depend  on  hk,  the  conduction  across  the 
feature  of  interest,  as  compared  with  Equation  (4) .  Without  the 
added  resistance  of  l/hk,  Equation  (8)  is  less  stable  than  Equation 
(4)  . 


Estimated  surface  temperatures  for  small  database  features  are 
obtained  from  analyzing  the  network  in  Figure  5b.  For  small 
features,  the  steady-state  surface  temperature  is  the  same  as  that 
in  Equation  (4)  with  the  exception  that  hk  =  d/k  instead  of  d/(2k). 
The  difference  in  the  2  representations  is  that  in  Equation  (4) 
small  objects  are  considered  to  be  attached  to  a  base  support  with 
temperature  T^^.  In  Figure  5b,  the  capacitance  network  for  small 
features,  an  extra  path  is  included  (from  Fig.  5a)  for  a  base  of 
support.  The  conductance  between  the  object  and  its  support  is  hk 
=  d/k. 

Figure  6  shows  a  plot  of  Equation  (8)  for  both  large  and  small 
features  across  a  daily  cycle  with  zero  wind  speed  and  rain  (i.e.. 
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Figure  5a.  Circuit  Representation  with  Capacitance 
for  Massive  Features. 


>.  Circuit  Representation  with  Capacitance 
for  Small  Features. 
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Figure  6.  Steady-State  Temperature  Estimates 
for  Capacitance  Model. 


both  hc  and  Rc  =  0) .  Note  the  large  variability  in  feature  surface 
temperature  across  time  for  large  features  (asphalt,  water,  and 
soil) .  This  variability  is  because  these  features  are  no  longer 
attached  to  a  base  support  as  in  Figure  2.  Surface  temperatures 
for  small  features  (i.e.,  the  bridge  crane  and  the  evergreen  tree) 
in  Figure  6  are  not  as  variable  across  a  diurnal  cycle  and,  in 
fact,  the  surface  temperature  of  the  bridge  crane  directly  mimics 
that  of  air  temperature  as  previously  found  (see  Fig.  1) . 

As  can  be  seen  by  comparing  Figure  6  with  Figure  1,  use  of  the 
capacitance  model  shown  in  Figure  5a  will  yield  predicted  feature 
surface  temperatures  for  large  features  with  a  much  wider  range 
than  the  original  circuit  in  Figure  2. 

The  complete  surface  temperature  solution  for  constant  radiant 
input  proposed  by  Hirsch  for  massive  bodies  consisted  of  a 
combination  of  a  steady-state  or  forced  component  and  a  time- 
varying  or  natural  component: 


T  = 


ss 


(T^-T 


- 1 
RC 


ss  ■‘oc; 


(1  cf/2 Jc  ^ 
R 


O 

Celsius 


(9) 


The  passage  of  time  in  seconds  is  denoted  as  t.  Tss  represents  the 
steady-state  temperature  component  described  previously  when  the 
capacitor  is  replaced  by  an  open  circuit.  Toc  represents  the 
temperature  of  the  database  feature  at  time  t  =  0  or  the  initial 
condition  of  the  feature.  The  resistance,  R,  denotes  the 
resistance  halfway  across  the  feature  itself,  R,  =  (d/2k) ,  in 
serial  with  the  remainder  of  the  resistance  throughout  the  circuit, 
R2  =  (l/(hc  +  Rc  +  4.62e))  f°r  large  features  or  R2  =  (l/(hc  +  Rc  + 
d/k  +  4.62  0)  for  small  objects  so  that  R  =  R1  +  R2 .  The  term 
d/(2k)/R  denotes  the  relative  weight  of  resistance  flowing  down  the 
path  of  the  capacitor  to  the  entire  resistance  in  the  circuit.  As 
this  weight  approaches  1  (i.e.,  there  is  infinitely  more  resistance 
down  the  path  of  the  capacitor  than  the  other  paths) ,  no  current 
will  flow  down  the  path  of  the  capacitor,  the  term  (1  -  d/(2k)/R) 
approaches  zero,  and  the  time-varying  component  in  Equation  (9) 
approaches  zero.  Capacitance  or  C  in  Equation  (9)  is  the  product 
of  the  density  and  specific  heat  of  the  material  composing  the 
feature  and  one-half  the  depth  of  the  feature.  Therefore,  thermal 
capacitance  is  a  linear  function  of  the  thickness  of  the  feature. 

Components  of  Equation  (9)  can  be  recombined  in  order  to  view 
feature  surface  temperature,  T,  as  a  linear  combination  of  the 
steady-state  temperature,  Tss,  and  the  initial  temperature,  T^. 
This  restructuring  is  as  follows: 

T  =  Tss  -  (Tss  -  T^wft)  =  (1  -  w(t))Tss  +  wft)^  (10) 
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where 


w ( t )  =  e"<t/<RC>>  ( l  -  d/(2k)/R) 


and  w(t)  represents  a  time-varying  weight  plotted  in  Figures  7  and 
8.  W(t)  decays  to  zero  for  all  features  as  time  increases.  This 
increase  forces  feature  surface  temperature  (T)  to  migrate  toward 
Tss  as  time  moves  forward  from  the  initial  condition. 

Individual  feature  material  properties  contributing  to  w(t) 
are  feature  emissivity,  thermal  conductivity,  material  density, 
heat,  and  material  depth  or  thickness.  For  the  sample  features, 
the  first  4  coefficients  are  provided  in  Table  1.  These 
coefficients  are  considered  constants;  however,  they  can  vary 
somewhat  as  a  function  of  temperature,  humidity,  etc.  The  final 
coefficient,  d  or  feature  thickness,  has  been  mentioned  previously 
as  a  source  of  concern  within  the  model,  affecting  resistance  and 
capacitance  in  the  model.  Feature  thickness  must  be  estimated 
(typically  by  the  database  modelers)  in  the  database  for  small 
features  and  through  computational  estimation  procedures  for 
massive  features  (Suits,  1985).  For  massive  features,  feature 
thickness  in  meters  is  estimated  as: 


86400klT 

nED 


(11) 


where  k  is  thermal  conductance,  E  is  specific  heat,  and  D  is  the 
material  density  for  the  database  feature  of  interest.  As 
mentioned  previously,  the  assumptions  accompanying  Equation  (11) 
may  be  somewhat  suspect  for  the  purpose  of  our  application.  The 
thickness,  d,  computed  in  Equation  (11)  for  massive  features, 
represents  the  maximum  thickness  of  material  involved  in  the 
diurnal  heat  transfer  such  that  the  temperature  at  this  depth 
decays  to  1/e  or  .368  of  the  surface  temperature  (see  Suits,  1985). 
For  the  5  sample  features,  estimates  of  d  are  provided  in  Table  1. 
The  computed  thickness  of  the  massive  features  (aircraft  parking, 
soil,  and  water)  is  greater  than  the  thickness  of  the  small 
features  (bridge  crane  and  evergreen  tree) .  It  is  not  clear, 
however,  whether  the  relative  differences  in  the  thicknesses  of  the 
materials  are  appropriate. 

Increments  in  the  thickness,  d,  of  a  feature  increases  (a)  the 
resistance  of  the  feature  to  thermal  flux,  and  (b)  the  capacitance 
of  the  feature.  These  increases,  in  turn,  have  an  effect  on  both 
the  steady-state  temperature,  Tss,  and  the  time-varying  component, 
(Tss  ~  Toe)  e-<t/(RC))  ( 1  -  d/(2k)/R). 

Increases  in  d  cause  a  decrease  in  the  conductance,  hk,  from 
Equation  (4)  ,  making  the  steady-state  temperature  of  a  feature  more 
dependent  on  input  from  the  sun,  sky,  rain,  and  air  as  opposed  to 
any  base  support  to  which  the  feature  is  attached. 
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(Long  Time  Scale) . 


Table  1.  Sample  Feature  Material  Coefficients 


Bridge 

crane 

Aircraft 

parking 

Soil 

Fresh 

water 

Evergreen 

tree 

Feature  Type 

small 

massive 

massive 

massive 

small 

Estimated 

Thickness (meters) 

.038 

— 

— 

— 

.006 

Computed 

Thickness (meters) 

.263 

.092 

.063 

— 

Density 

7800 

500 

466 

1000 

500 

Thermal 

Conductivity  (k) 

40 

1.047 

.1260 

.6 

.15 

Specific  Heat  (E) 

470 

830 

879 

4180 

2400 

Capacitance  (C) 

69654 

54572 

18842 

131670 

3600 

Feature  Conductance 

2k/d 

2105 

8 

2.7 

19 

50 

Capacitance  Circuit 

Resistance  (R) 

00142 

.355 

.6 

.28 

.054 

Massive  objects 

i  do  not 

resist  the 

input  flux  more  ; 

because  of 

their  thickness  or 

size  as 

i  the  model 

from 

Equation 

(9)  would 

currently  indicate.  These  objects  distribute  the  thermal  energy 
throughout  the  material.  This  notion,  however,  is  indicative  of 
greater  capacitance  and  should  be  reflected  through  the  capacitance 
parameter.  Soil,  for  example,  is  a  massive  object  which,  although 
it  does  not  conduct  thermal  energy  well,  has  great  capacity  for 
distributing  the  thermal  energy  it  does  conduct.  Six  to  twelve 
feet  below  the  surface,  the  temperature  of  the  earth  remains 
constant  at  approximately  60  °F.  The  temperature  gradient  from  the 
surface  of  the  earth  to  a  depth  of  6  feet  is  continuous.  This 
continuity  constraint  limits  the  effect  of  radiant  energy  or 
atmospheric  air  temperature  in  modifying  the  temperature  at  the 
soil's  surface.  In  the  thermal  modeling  of  soil,  this  constraint 
could  be  modeled  by  placing  a  resistor  of  value  d/k  where  d  =  2 
meters  and  k  is  .  12  60  and  a  voltage  source  of  60  °F  or  15.6  °C  below 
the  capacitor  representing  soil  in  the  circuit. 
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The  thickness  of  the  soil  layer,  however,  should  not  affect 
soil's  ability  to  receive  radiant  flux.  Soil,  however,  also  has 
low  thermal  conductance,  k,  as  opposed  to  metal.  Thermal 
conductance  is  inversely  related  to  the  resistance  of  a  material  to 
thermal  flux  as  one  may  attest  to  by  comparing  the  feel  of  metal 
with  that  of  a  piece  of  soil  on  a  hot  summer  day.  The  use  and 
placement  of  an  object's  resistance,  d/k,  in  our  circuit  model  is 
misleading.  Placing  this  resistance  in  front  of  our  capacitor 
prevents  thermal  flux  from  entering  the  capacitor.  As  can  be 
inferred  from  Fourier's  Law,  this  resistance  should  actually  apply 
only  to  objects  on  the  opposite  side  of  the  feature  of  interest. 
This  resistance  is  one  of  the  problems  of  representing  thermal 
gradients,  which  are  continuous  in  space,  with  discrete  circuit 
models.  This  source  of  error  plays  a  major  role  in  the  thermal 
model . 

To  estimate  the  initial  temperature,  Toc,  from  Equation  (9) , 
Hirsch  suggested  that  the  software  could  be  run  at  20-min  intervals 
1  day  prior  to  the  estimate  of  Toc  actually  desired.  T^  should  be 
set  to  zero  for  the  initial  run  ("i"  =  1)  with  subsequent  values  of 
Toc  on  computer  run  "i  +  l"  set  to  the  final  estimate  of  T  for  the 
respective  feature  from  run  "i." 

In  Equation  (9)  ,  the  final  feature  surface  temperature,  T, 
will  be  the  steady-state  temperature  Tss  modified  by  a  time-varying 
component  which  will  be  a  percentage  of  the  difference  between  Tss 
and  JL-  The  decaying  exponential  in  the  time-varying  component, 
e  (t  (where  t  is  expressed  in  seconds) ,  places  less  weight  on  the 
time-varying  component,  (Tss-Toc)  ,  as  time  passes.  The  other  weight 
in  the  time-varying  component,  (1  -  d/2k/R) ,  denotes  relative 
resistance  down  the  path  of  the  feature  of  interest  to  the  entire 
circuit  and  varies  between  zero  and  one.  If  the  resistance  d/(2k) 
down  the  path  of  the  capacitance  is  quite  large  relative  to  the 
total  resistance  of  the  circuit,  very  little  flux  will  flow  into 
the  capacitor  and  there  will  be  no  charging  of  the  capacitor,  or 
equivalently,  no  thermal  heating  of  the  feature.  In  this  instance, 
the  term  (1  -  d/(2k)/R)  will  be  approximately  zero  and  the  feature 
surface  temperature  will  always  be  the  steady-state  temperature  Tss. 

The  thermal  heat  transfer  process  is  continuous  in  time  and 
three-dimensional  space.  The  electrical  circuit  is  composed  of 
discrete  elements.  For  example,  in  the  circuit  model  the 
resistance  to  flux  input  at  the  feature  capacitor  is  d/(2k).  This 
simple  representation  means  that  the  flux  travels  through  half  of 
the  material,  d/2,  before  it  reaches  the  capacitor.  As  the  depth 
of  the  material  d  increases,  less  flux  reaches  the  capacitor.  This 
is  a  discrete  representation  of  a  continuous  process.  In  an  actual 
feature  of  reasonable  size,  a  temperature  gradient  will  exist 
across  the  material  of  interest  and  capacitance  exists  throughout 
the  material,  not  just  at  a  depth  of  d/2.  For  example,  large  metal 
objects  heated  on  one  side  may  exhibit  a  large  temperature  gradient 
from  the  heated  or  sunny  side  to  the  underside.  Soil,  which  is  one 


19 


of  the  sample  features  of  interest  in  this  report,  is  a  poor 
conductor  of  thermal  energy  and  will  not  heat  up  like  metal  in 
response  to  the  sun.  Six  feet  below  the  surface,  soil  temperature 
remains  stable  at  approximately  60  °F.  To  properly  model  surface 
temperature,  this  gradient  from  below  the  surface  to  the  surface  of 
the  soil  should  be  properly  modeled.  Although  the  example 
described  here  can  be  more  accurately  modeled  through  discrete 
circuitry  than  is  presently  the  case,  it  would  be  necessary  to 
implement  such  improvements  on  a  feature-by-feature  basis. 

The  characteristics  of  the  time-varying  weights,  e  and 

(1  -  d/(2k)/R),  may  be  analytically  studied  to  determine  their 
general  influence  on  feature  surface  temperature  estimation,  or  for 
approximating  T0(;,  as  Hirsch  suggested.  If  either  of  the  weights 
or  the  combination  of  weights  approaches  zero,  the  time-varying 
component  of  the  surface  temperature  equation  will  not  contribute 
to  Equation  (9)  .  Figures  7  and  8  show  the  multiplicative 

combination  of  weights  (i.e.,  e  (t/c  (1  -  d/(2k)/R))  for  the  time- 
varying  components  of  feature  surface  temperatures  for  the  5  sample 
features  (bridge  crane,  aircraft  parking  area,  soil,  water,  and  an 
evergreen  tree) .  Note  that  the  aircraft  parking  area,  soil,  and 
water  are  considered  massive  objects  while  the  bridge  crane  and 
evergreen  tree  are  small  objects. 

At  ^ime  zero  in  Figures  7  or  8,  the  decaying  exponential  is 
simply  e  =1  and  the  multiplicative  weight  for  each  sample  feature 
becomes  (1  -  d/(2k)/R) .  According  to  the  model,  water  provides  the 
least  amount  of  resistance  to  heat  with  a  weight  in  the  neighbor¬ 
hood  of  .76  at  a  time  of  zero.  This  finding  is  shown  in  Figures  7 
and  8.  As  mentioned  previously,  the  component  (l-d/(2k)/R) 
represents  the  resistance  the  radiant  flux  meets  from  the  feature 
of  interest  relative  to  the  components  in  the  circuit.  Large 
features  might  show  more  capacitance  but  not  necessarily  more 
resistance  to  flux  at  the  surface.  For  example,  in  Figure  5,  water 
shows  more  resistance  to  incoming  flux  than  does  the  evergreen 
tree.  This  is  not  logical  as  we  would  expect  that  water  conducts 
thermal  radiation  more  readily  than  an  evergreen  tree. 

The  decay  of  the  curves  in  Figures  7  and  8  is  determined  by 
their  time  constants,  1/(RC),  from  the  exponent  in  Equation  (9). 
The  larger  the  time  constant,  the  more  rapid  is  the  decay  of  the 
curve  in  Figures  7  and  8.  Smaller  time  constants  (larger  RC 
values)  represent  capacitors  which  take  longer  to  store  up  energy 
as  well  as  dissipate  energy.  Thus,  we  might  expect  that  massive 
objects  would  have  smaller  time  constants  (larger  RC  values) 
because  thermal  energy  should  build  up  as  well  as  dissipate  more 
slowly  over  time  for  these  objects.  This  premise  holds  true  as 
expressed  in  Figures  7  and  8.  In  Figure  8,  decay  rates  of  the 
sample  massive  features  are  available.  Soil  has  the  fastest  decay, 
followed  by  the  aircraft  parking  area  and  water.  These  storage  and 
dissipation  rates  do  not  meet  our  intuitive  expectations. 
Logically,  we  would  expect  soil  to  have  the  slowest  uptake  and 
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release  of  heat,  followed  by  water  and  then  asphalt  from  the 
aircraft  parking  area.  Mathematically,  these  expectations  would  be 
reflected  in  the  RC  combinations  with: 

R^water  <  R^aircraft  p  <  R^soit  (12) 

or  equivalently 

1/RC«ater  >  i/Rc  aircraft  p  >  l/RCsoil.  (13) 

In  the  time  constants  1/ (RC)  we  can  compare  both  the  values  of 
R  and  C  for  small  versus  large  features.  As  expressed  previously, 
R  =  d/(2k)  +  l/(4.6e  +  hc  +  Rc)  for  large  features  and  R  =  d/(2k) 
+  1/(4.626  +  hc  +  Rc  +  k/d)  for  small  features  and  only  d  and  k,  the 
material  depth  and  conductivity,  change  significantly  across 
features.  When  there  is  no  rain  or  wind,  l/(4.6e  +  hc  +  Rc)  is 
approximately  .23  and  R  becomes  approximately  d/(2k)  +  .23  for 

large  features.  For  small  features,  R  is  approximately 
[(3k  +  4 . 37d)  /  (k  +  4 . 37d)  ]  X  d/(2k).  Estimates  of  R  under 

conditions  of  no  rain  or  wind  are  given  in  Table  1  for  the  sample 
features.  Capacitance  or  C  in  the  term  1/ (RC)  is  the  product  of 
material  density,  the  specific  heat  of  the  material,  and  one-half 
the  depth  of  the  material.  Capacitance  for  each  of  the  sample 
features  is  given  in  Table  l.  As  can  be  seen  from  Table  1,  soil 
has  the  smallest  estimated  capacitance.  This  finding  is  not  as 
expected  and,  in  fact,  the  ground  or  soil  stores  and  releases 
thermal  energy  much  more  slowly  than  small  bodies  of  water.  The 
value  of  d,  the  thickness  of  the  soil,  could  be  increased  to 
increase  C,  the  capacity  for  soil,  relative  to  the  other  sample 
features.  Estimation  of  material  ^  ._ckness  for  massive  features 
was  accomplished  through  Equaticn  (11),  as  previously  mentioned, 
and  these  estimates  are  provided  in  Table  1.  In  the  present  model, 
capacitance  is  a  linear  function  of  the  depth  of  the  features  so 
that  doubling  the  r’epth  of  a  feature  would  also  double  the 
capacitance.  In  Table  l  we  see  that  the  depth  of  the  massive 
features  (soil,  aircraft  parking  area,  and  water)  is  not 
necessarily  greater  than  the  depth  of  the  small  features  (bridge 
crane  and  evergreen  tree) . 

Intuitively,  then,  the  relative  ordering  of  the  thicknesses 
for  the  5  features  in  Table  1  is  questionable.  Because  massive 
feature  depths  are  estimated  from  Equation  (11) ,  the  use  of 
Equation  (11)  in  estimating  material  depth  for  use  in  capacitance 
is  considered  questionable. 

Given  this  consideration  of  estimating  feature  thickness,  we 
can  return  to  the  problem  of  approximating  Equation  (9)  through 
numerical  iteration.  Hirsch  suggested  that  feature  temperature 
estimates,  T,  be  computed  every  twenty  minutes  with  TK  initially 
set  to  zero.  On  each  subsequent  iteration,  then,  T^  is  set  to  the 
value  of  T  computed  on  the  previous  iteration.  From  Figure  7  we 
can  estimate  that,  at  one-third  of  an  hour,  the  time-varying 
weights  for  the  small  features,  the  bridge  crane  and  the  evergreen 
tree,  are  approximately  zero.  This  results  in  the  estimate  of  T 
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always  being  approximately  Tss.  For  the  massive  features  of  water, 
the  aircraft  parking  area  and  soil,  the  weights  at  t  =  one-third 
hour  (from  Fig.  8)  are  approximately  .76,  .58,  and  .34, 
respectively.  New  estimates  of  soil  surface  temperature,  then,  are 
66%  of  the  current  feature  steady-state  temperature  plus  34%  of  the 
surface  temperature  of  the  soil  from  20  min  prior. 

We  can  ask  how  this  reformulation,  using  the  time-varying 
component,  modifies  feature  surface  temperatures  across  a  diurnal 
cycle.  Our  general  goal  (looking  back  to  Fig.  1)  is  to  soften  the 
extreme  temperature  variation  for  features  such  as  soil.  The 
iterative  approach  suggested  by  Hirsch  (iteration  every  20  minutes 
to  update  feature  surface  temperature)  leads  to  the  equation: 

j  -  0 


where  T,  is  the  feature  surface  temperature  on  iteration  "i,"  Ti+W 
is  the  feature  surf a 
is  the  steady-state 
denotes  the  weight, 

(10) .  Equation  (14)  shows  that  the  feature  surface  temperature  at 
any  time  is  a  weighted  combination  of:  (a)  the  feature  surface 
temperature  at  any  previous  time  plus,  and  (b)  a  sum  of  steady- 
state  temperatures  between  the  initial  period  and  the  present  time. 
Figures  9  and  10  show  estimates  of  Equation  (11)  for  soil  (P  =  w(t) 
*  .34)  and  water  (P  =  w(t)  «  .76),  respectively.  The  steady-state 
temperature  in  Figures  9  and  10  is  a  sine  wave  used  to  approximate 
diurnal  variation.  In  both  Figures  9  and  10,  the  initial  feature 
surface  temperature  at  time  zero  or  "iw  =  0  is  X,  =  0,  as  Hirsch 
suggested.  Iterative  updates  in  Equation  (14)  were  made  every  20 
min.  As  can  be  seen  from  Figures  9  and  10,  the  iterative  estimates 
quickly  fall  into  line  with  the  sine  wave  as  time  passes.  As  the 
value  of  P  increases  from  .34  (soil)  (Fig.  9)  to  .76  (water)  (Fig. 
10) ,  the  iterative  estimate  tends  to  lag  more  behind  the  steady- 
state  input  but  the  effect  is  still  minimal. 

As  mentioned  previously,  our  goal  in  adding  heat  storage 
capacity  to  the  model  was  twofold.  First,  it  was  necessary  to 
soften  extreme  surface  temperature  values  across  the  daily  heating 
and  cooling  cycle,  as  exemplified  by  soil  in  Figure  1.  Second,  it 
was  desirable  to  create  a  time  lag  such  that  massive  features  heat 
slowly  over  the  day  and  release  heat  slowly  as  the  air  temperature 
cools  after  the  sun  goes  down.  The  surface  temperatures  predicted 
by  the  present  iterative  model  will  not  perform  either  of  these 
functions  because  the  iterative  estimate  closes  in  on  the  steady- 
state  temperature  too  quickly  in  time. 


.«=  temperature  on  iteration  "l  +  N,M  Tss(1-  ♦  j  ♦  1} 
temperature  on  iteration  "i  +  j  +  l,"  and  P 
w(t)  =  e'(  ”(1  -  d/(2k)/R),  from  Equation 
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Figure  9.  Iterative  Capacitance  Model  Predictions 
for  Weight  =  . 34 . 


CIRCUIT  SOLUTION  FOR  TRANSIENT  INPUT 

The  surface  temperature  solution  provided  by  Hirsc^/(R£ijom 
Equation  (9)  in  the  previous  section  •  T  Tss  “  ( Tss  -  Tpc )  e 
(1  -  d/(2k)/R),  is  a  solution  to  a  first-order  differential 

equation  with  constant  input  flux.  Flux  input  from  the  sun  and  sky 
does  vary  in  a  sinusoidal  fashion  across  a  diurnal  cycle;  however, 
the  solution  provided  above  is  an  approximation  for  slow  changes  in 
input  flux  relative  to  the  response  of  the  storage  element (s).  In 
this  section,  a  more  complete  solution  is  provided  to  the  network 
from  Figure  2  in  order  to  compare  performance  of  the  more  complete 
solution  with  the  approximation  from  the  previous  section. 

Flux  input  from  the  sky  and  sun  minus  the  constant  (316e)  from 
the  linear  black  body  approximation  may  be  represented  by  a 
sinusoidal  waveform  of  the  form: 

Ps„  +  Psky  -  316e  =  Asin (wt  -  0)  +  B  (15) 

where  t  is  expressed  in  seconds  and  t  =  0  designates  midnight,  w  = 
2rrf  -  2ti  (1  cycle/day)  =  2ir  (1  cycle/86,400  seconds)  =  7.27  X  10 
radians/second,  0  =  n/2  designating  2  a.m.  as  the  time  of  lowest 
flux  input  from  the  sun  and  sky,  and  B  -  A  being  the  minimum  flux 
input  (occurring  at  2  a.m.)  while  A  +  B  denotes  the  maximum  flux 
input  (occurring  at  2  p.m.)  across  the  diurnal  cycle.  Figure  11 
shows  a  sample  plot  of  Equation  (15). 

Employing  Equation  (15)  from  above  in  the  circuit  in  Figures 
5a  and  5b  yields  a  nonhomogeneous  first-order  differential 
equation.  Using  circuit  theory  properties  of  linearity  and 
superposition,  the  solution  for  T,  the  feature  surface  temperature, 
may  be  arrived  at  by  solving  the  circuit  for  the  response  to  the 
sinusoidal  input  (Asin  (wt  +  0))  and  the  constant  input  (B) 
separately,  and  then  adding  the  two  responses  together.  The 
feature  surface  temperature  solution  for  the  constant  input  is  the 
solution  to  Equation  (9) ,  the  homogeneous  first-order  differential 
equation  provided  in  the  last  section  of  the  form: 

T1  =  Tss  -  (Tss  -  T0c)e'<t/<RC>>(1  -  d/(2k)/R)  °Celsius  (16) 

All  of  the  parameters  in  Equation  (16)  are  the  same  as  in  Equation 
(9)  with  the  exception  of  the  steady-state  temperature  which  now 
takes  the  value: 


Tss 

for  massive 


B  +  hcTajr  +  RcTrajn 
4 . 6  €  +  hc  +  Rc 
features,  and 


“Celsius  (17a) 
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iance  Estimates 


4 . 6  e  +  hc  +  Rc  +  (k/d) 


“Celsius  (17b) 


for  small  features.  All  parameters  in  Equations  (17a)  and  (17b) 
are  identical  to  Equation  (8)  with  the  exception  of  B,  the  constant 
flux  input  from  Equation  (15)  . 

The  response  to  the  sinusoidal  input,  Asin(wt  +  0),  may  be 
obtained  using  the  Laplace  Transform,  yielding  a  feature  surface 
temperature 

T2  =  Rsin(wt  +  0+0,)  (18) 


where 


R= - 

(4.62  e+hc+Rc2) 
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for  massive  features,  and 


R= 


^^c^'aii^^-c'^'zain*'  ^  ^baae 
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for  small  features.  T2  is  a  sinusoidal  waveform  which  follows  the 
sinusoidal  form  of  the  input  flux  from  the  sun  and  sky  with  the 
exception  of  a  change  in  magnitude  and  a  phase  lag,  denoted  by  R 
and  0,,  respectively. 

Figure  12  shows  a  plot  of  Equation  (18)  for  the  5  sample 
features.  Note  that  the  temperatures  in  Figure  12  fluctuate  around 
0  °F.  The  DC  (direct  current)  or  average  temperature  component  is 
not  included  in  Equation  (18)  nor  in  Figure  12.  Figure  12  shows 
only  the  variability  in  temperature  across  a  diurnal  cycle. 
Certainly  the  range  in  temperatures  across  a  single  day  for 
individual  features  in  Figure  12  is  too  large.  For  example,  in 
Figure  12  evergreen  trees  have  a  range  of  approximately  120  °F  and 
water  has  a  range  of  approximately  60  °F.  This  shortcoming 
directly  follows  from  the  resistive  inputs  to  the  magnitude,  R,  in 
Equation  (18)  and  can  be  attributed  to  inadequacies  in  the  circuit 
model.  The  conclusion  here  is  that  another  parallel  conduit  or 
path  in  the  circuit  model  is  required  to  carry  more  thermal  flux 
away  from  the  feature  of  interest. 

The  important  parameters  in  Figure  12  are  the  magnitudes  of 
the  sine  waves,  R,  which  represent  minimum  and  maximum  daily  heat 
values,  and  the  phase  shift,  0V  which  represents  time  lag  in  the 
buildup  and  release  of  heat  throughout  a  daily  cycle.  For  the  5 
sample  features,  these  values  are  presented  in  Table  2. 

Estimates  of  magnitudes  in  Table  2  are  much  like  estimates  of 
Tss  in  the  initial  part  of  this  report.  If  we  include  a  base  of 
support  in  the  circuit  model  as  we  have  done  with  small  features, 
the  magnitudes  will  be  dependent  upon  the  base  support  temperature 
to  the  extent  that  feature  conductance,  hk  =  k/(2d),  outweighs  all 
other  conductances.  For  the  bridge  crane,  hk  is  quite  large  (see 
Table  1) ,  forcing  the  magnitude,  R,  to  be  approximately  the  base 
temperature  (air  temperature) .  For  the  evergreen  tree,  hk  is  much 
more  in  line  with  other  conductance  values  in  the  circuit,  and  R 
tends  to  be  an  averaging  of  air  temperature  with  irradiance  from 
the  sun  and  sky.  For  the  massive  features,  the  estimated 
temperature  magnitudes  are  dependent  upon  the  irradiance  from  the 
sun  and  sky. 
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Figure  12.  Phasor  Solution  for  Capacitance  Model. 


Table  2.  Sample  Feature  Magnitude  and  Phasor  Coefficients 


Bridge 

crane 

Aircraft 

Darkina 

Soil 

Water 

Evergreen 

tree 

Magnitude  (R) 
“Fahrenheit 

air 

temp. 

82.1 

81.4 

76.1 

70  °F* 

Phase  (0,) 

.27° 

19.7° 

9.2° 

29.3° 

.48° 

"Approximately  the  average  of  large  feature  magnitudes  with 
the  value  of  air  temperature — in  this  instance  70  °F. 


For  both  small  and  massive  features,  the  estimated  magnitudes 
in  Table  2  are  much  too  large.  Recall  that  the  results  from  Figure 
12  and  Table  2  indicate  only  the  modulation  of  daily  feature 
surface  temperature  around  an  average  level  which  is  denoted  by 
Equation  (16) .  Feature  surface  temperatures  will  not  vary 
throughout  a  daily  cycle  by  100  °F  (-20  °F  from  the  average 
temperature  at  2  a.m.  and  +80  °F  at  2  p.m.).  This  shortcoming  in 
the  circuit  model  has  pervaded  all  of  our  analyses.  Only  by 
including  more  resistance  in  the  circuit  or  decreasing  the 
modulation  in  the  daily  irradiance  from  the  sun  and  sky  can  this 
problem  be  rectified. 

Phase  angles  or  the  lag  in  heat  accumulation  by  the  features 
is  also  problematic,  as  shown  in  Figure  12  and  Table  2.  Table  2 
shows  that  water  heats  and  cools  more  slowly,  followed  by  the 
aircraft  parking  area,  soil,  and  finally  the  2  small  features. 
These  findings  are,  as  discussed  in  the  previous  section,  directly 
related  to  the  estimated  capacitance  for  these  features  as  shown  in 
Table  1.  For  small  features,  the  lag  is  approximately  zero  because 
their  conductances,  hk  (see  Table  1)  ,  are  large  relative  to  all 
other  factors  in  the  denominator  of  the  phase  angle  calculation 
(see  Equation  (18) ) .  Of  course,  these  conductance  values  are 
linearly  related  to  d,  the  estimated  feature  thickness  (i.e.,  hk  = 
2k/d) ,  and  so  would  be  quite  sensitive  to  any  changes  in  estimated 
feature  thickness. 

The  revised  capacity  model  in  this  section  provided  a  better 
mechanism  for  allowing  heat  to  be  stored  and  released  from 
individual  features  over  time.  Individual  feature  material 
parameters  (e.g.,  feature  thickness)  and  the  circuit  representa¬ 
tion,  however,  are  still  problematic,  preventing  the  model  from 
producing  reasonable  surface  temperatures  across  a  daily  cycle. 
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CONCLUSIONS 


Heat  transfer  between  media  is  a  complex  process.  The 
idiosyncratic  nature  of  the  process  makes  it  difficult  to  develop 
a  general  model  for  predicting  feature  surface  temperatures  as  a 
function  of  environmental,  geographic,  temporal,  and  feature 
material  input. 

In  the  current  IR  predictive  model  in  use  at  AL/HRA, 
shortcomings  in  a  steady-state  thermal  predictive  model  led  to 
large  inaccuracies  in  predicting  feature  surface  temperature  for 
database  features  such  as  soil.  Factors  such  as  type  of  material, 
feature  thickness,  the  attachment  of  a  feature  to  a  base  of 
support,  and  the  extent  of  the  support  are  critical  to  predicting 
heat  transfer  and,  consequently,  feature  surface  temperature. 
Representing  the  continuous  process  of  heat  transfer  in  three- 
dimensional  space  by  resistors,  capacitors,  and  current  and  voltage 
sources  in  an  electrical  circuit  necessarily  involves  simplifi¬ 
cation  of  the  heat  transfer  process  and  the  model  suffers  in  the 
simplification  process  as  discussed  throughout  this  report. 

Inclusion  of  heat  capacity  or  the  ability  to  store  heat  had  no 
significant  effect  on  softening  extreme  diurnal  feature  surface 
temperatures  or  providing  a  temperature  lag  throughout  a  diurnal 
cycle.  Of  course  these  findings  are  dependent  upon  the 
mathematical  repi _  -ntation  of  capacitance  and  simply  infer  that 
the  model  is  ir  z^c-guate.  Development  of  a  more  valid  capacitance 
structure,  however,  would  have  to  be  couched  iu  the  circuit 
representation  and,  from  the  present  analyses,  it  is  evident  that 
a  more  complex  circuit  representation  is  required. 

The  extreme  feature  surface  temperatures  across  a  daily  cycle 
found  in  this  report  and  previous  reports  may  be  alleviated  by 
either — (a)  adding  more  resistive  elements  in  parallel  with  the 
current  circuitry  of  the  model  which  help  to  drain  irradiance  from 
the  sun  and  sky,  or  (b)  filtering  the  irradiance  from  the  sun  and 
sky  to  decrease  the  total  amount  of  input  flux  to  the  circuit. 

The  relative  heating  and  cooling  of  individual  features  across 
a  daily  cycle  may  be  represented  by  including  capacitors  in  the 
circuit  model.  However,  given  these  more  complex  circuit 
representations,  the  predictions  from  the  circuit  representation 
are  still  limited  by  the  accuracy  of  individual  material 
coefficients.  For  example,  relative  lag  time  in  the  heating  and 
cooling  of  the  sample  features  in  this  report  are  entirely 
predicted  by  the  electrical  circuit  capacity  parameters.  The 
orderings  of  these  individual  capacitance  parameters  for  sample 
features  in  this  report,  however,  were  not  as  intuitively  expected 
(e.g.,  we  would  intuitively  expect  water  to  heat  and  cool  more 
quickly  than  soil  but  their  estimated  capacitance  parameters  from 
this  report  contradict  this  expectation) . 
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